Investigating the shared genetic architecture between primary sclerosing cholangitis and inflammatory bowel diseases: a Mendelian randomization study

Background Several studies have found that primary sclerosing cholangitis (PSC) and inflammatory bowel disease (IBD) are closely associated. However, the direction and causality of their interactions remain unclear. Thus, this study employs Mendelian Randomization to explore whether there are causal associations of genetically predicted PSC with IBD. Methods Genetic variants associated with the genome-wide association study (GWAS) of PSC were used as instrumental variables. The statistics for IBD, including ulcerative colitis (UC), and Crohn’s disease (CD) were derived from GWAS. Then, five methods were used to estimate the effects of genetically predicted PSC on IBD, including MR Egger, Weighted median (WM), Inverse variance weighted (IVW), Simple mode, and Weighted mode. Last, we also evaluated the pleiotropic effects, heterogeneity, and a leave-one-out sensitivity analysis that drives causal associations to confirm the validity of the analysis. Results Genetically predicted PSC was significantly associated with an increased risk of UC, according to the study (odds ratio [OR] IVW= 1.0014, P<0.05). However, none of the MR methods found significant causal evidence of genetically predicted PSC in CD (All P>0.05). The sensitivity analysis results showed that the causal effect estimations of genetically predicted PSC on IBD were robust, and there was no horizontal pleiotropy or statistical heterogeneity. Conclusions Our study corroborated a causal association between genetically predicted PSC and UC but did not between genetically predicted PSC and CD. Then, we identification of shared SNPs for PSC and UC, including rs3184504, rs9858213, rs725613, rs10909839, and rs4147359. More animal experiments and clinical observational studies are required to further clarify the underlying mechanisms of PSC and IBD. Supplementary Information The online version contains supplementary material available at 10.1186/s12876-024-03162-6.


Introduction
Inflammatory bowel disease (IBD) is a chronic intestinal disorder with unknown etiology.Many studies point to the presence of genetic predisposition, intestinal mucosal immune system dysfunction, and microbiota imbalance [1] in the occurrence and progression of IBD.Ulcerative colitis (UC) and Crohn's disease (CD) are two main typical subtypes of IBD.The incidence of IBD has risen over the past decade in Asia.Predictably, the prevalence of IBD will significantly in the future, following an aging population [2].
Patients with IBD not only suffer a significant reduction in their quality of life but also causes substantial costs in health care due to its high prevalence [3].Chronic IBD is restricted to the gut, but also in the extraintestinal organs in many patients [4,5].This phenomenon is called extraintestinal manifestations (EIM) of IBD.EIM frequently affects joints [6], skin [7], eyes [8], lungs [9], pancreas [10], and liver [11].Primary sclerosing cholangitis (PSC) is important EIM in IBD patients [4].In clinical, about 70% of PSC patients are found to have underlying IBD [12][13][14].Genetic risk factors, environmental factors, activation of the immune system, and microbiota have been assumed that the factors relevant to the pathogenesis of EIMs [15,16].For PSC, the association with the activity of the underlying IBD is unclear [17].
PSC is a type of autoimmune liver disease characterized by multi-focal bile duct strictures and progressive liver disease [18].The prognosis of PSC was not satisfactory.Most patients ultimately require liver transplantation, after which disease recurrence may occur.However, without liver transplantation, the median survival time for PSC patients is 10 to 12 years [19].
Similar to IBD, the pathogenesis of PSC is also not well clear.However, the characteristic that PSC is often accompanied by IBD suggests that there may be a shared pathogenic gene or pathway between the PSC and IBD.Mendelian randomization (MR) renders us a novel way to study the connection between these two diseases.MR is a genetic epidemiological method, this method follows the Mendelian genetic law of "parental alleles are randomly assigned to offspring" [20].This method uses singlenucleotide polymorphisms (SNPs) as instrumental variables (IVs) to infer the potential causality of exposure and outcome.It's beneficial to minimize bias caused by confounding factors and reverse causality [21].Based on this, the MR method has been widely used to assess the causal relationship between traits and diseases or between diseases [22][23][24].To the best of our knowledge, there are no MR studies inferring the potential causal relationship of PSC with IBD to date.Therefore, we applied the MR method to examine whether the genetically predicted PSC is associated with IBD.

Study design
The overall design of Mendelian randomization analyses in this study is shown in Fig. 1.Briefly, (1) The selected instrumental variables were linked to exposure; (2) there were no inherent interactions between instrumental variables and confounder factors; (3) exposure was the only way by which instrumental variables can affect outcomes.The PSC served as the exposure, and UC served as the outcome.Since all datasets used in this study were based on public databases, no additional ethical approval was required.

GWAS data for PSC, UC, and CD
We gathered the summary statistics of PSC, UC, and CD, from the IEU Open GWAS project (https:// gwas.mrcieu.ac.uk/), all the cases there were defined on the basis of the International Classification of Diseases (ICD), and fulfilled the clinical diagnosis criteria for IBD and PSC.To be more specific, the sample sizes of datasets for PSC, UC, and CD, are 14,890 cases, 463,010 cases, and 212,356 cases, respectively.PSC has 7, 891, 603 SNPs, UC  1.

Instrumental variable selection
All statistical analyses were performed by the R packages: TwoSampleMR.First, we selected SNPs related to PSC at the genome-wide significance threshold with p< 5 × 10 -8 .Because strong linkage disequilibrium could lead to biased results.Second, the independence among the selected SNPs was evaluated according to the pairwiselinkage disequilibrium (r2 < 0.001, clumping window of 10,000 kb).When F-statistics were greater than 10, SNPs were considered powerful enough to mitigate the influence of potential bias.Third, we selected SNPs with F statistic >10 as IVs.

Statistical analysis
Based on the IVs, we performed an MR analysis to investigate the relationship between PSC and IBD.Five popular MR methods were used to analyze our data: MR Egger, Weighted median (WM), Inverse variance weighted (IVW), Simple mode, and Weighted mode.The IVW method is reported to be slightly more powerful than the others under certain conditions.
Cochran's Q statistics were used to perform heterogeneity, and p > 0.05 indicated no heterogeneity.Moreover, the MR-Egger method was used to determine the horizontal pleiotropy, MR-Egger at a p-value < 0.05 can imply the presence of horizontal pleiotropy.

Selection of instrumental variables
After a series of quality control steps as mentioned above, 18 SNPs were selected as IVs (Table 2).

Causality relationship between PSC and IBD
Among the five MR methods, the causal effects of genetically predicted PSC on UC and CD were inconsistent.The results of the MR analyses were shown in Table 3, genetically predicted PSC was positively associated with a risk of UC in our study, with a p-value of IVW method less than 0.05.However, we found no evidence supporting a causal association between PSC and CD.Previous research indicated that the genome-wide genetic correlation between PSC and UC was significantly greater than that between PSC and CD [25], similar to our results.
The scatter plots were used to show the single SNP effect and the combined effects of each MR method (Fig. 2).Forest plots and funnel plots of the causal effect are shown in Supplementary Figure 1.

Sensitivity analysis
We performed a leave-one-out sensitivity analysis, heterogeneity, and horizontal pleiotropy to further verify the reliability of our results.The results of sensitivity analysis showed that the causal effect estimation of this study was robust.The MR-Egger (Q p-value 0.137) and IVW methods (Q p-value 0.214) showed no statistical heterogeneity.Furthermore, no statistical horizontal pleiotropy was found in the horizontal pleiotropy of MR-Egger methods (P=0.719).The results of the sensitivity analysis are shown in Supplementary Figure 2.

Discussion
The etiology of PSC and IBD remains unclear, and there is a lack of effective treatment methods.Now, the main treatment methods for PSC include bile composition  modulators, immune modulators, anti-fibrotic, and regulation of the microbiome.However, further research is needed to determine whether these methods can delay its progression or improve transplant-free survival [26].
The same applies to the treatment of IBD.Although some new methods such as fecal transplantation, and small molecule drugs, applied to the treatment of IBD, satisfactory results have not been achieved in clinical yet [27][28][29].Therefore, it is crucial to investigate the relationship between PSC and the subtype of IBD.
Previous studies have suggested an association between IBD and PSC or PBC [30,31].PSC is a prototypic gut-liver axis disease.In the patients of PSC, gut microbiota could disrupt the intestinal barrier, leading to bacterial translocation and Th17 cell-driven liver damage [32].In contrast, the bile acid metabolizing enzyme CYP8B1 inhibits self-renewal of crypt based   intestinal stem cells through the accumulation of its product bile acid, hinders intestinal epithelial barrier repair, and exacerbates inflammatory response [33].These studies indicated a close correlation between intestinal diseases and liver diseases.As mentioned earlier, genetic predisposition plays a role in the occurrence and progression of IBD and PSC.The formation of serum antibodies is a way in which genetic factors affect the immune system.Multiple antibodies such as anti-Saccharomyces cerevisiae antibodies (ASCA), anti-neutrophil cytoplasm antibodies (ANCA) were upregulated in both autoimmune live diseases and IBD [34][35][36], and those antibodies may predict development of disease.The expression of common antibodies can also indicate a close relationship between the two diseases [37].
In this study, we used GWAS data to investigate the possible causal relationship and specific SNPs between PSC and IBD susceptibility, offering novel insights into the prevention and treatment of PSC and IBD.Multiple MR methods were employed to investigate the relationship between PSC and UC or CD, respectively.Four MR methods (Weighted median, Inverse variance weighted, Simple mode, and Weighted mode) indicated a significant relationship between PSC and UC.However, as for CD, there was no significant relationship between PSC and CD.Thus, we conclude that PSC has a significant relationship with UC but not CD.According this analysis, we also found the specific SNPs that are shared for PSC and UC (rs3184504, rs9858213, rs725613, rs10909839, and rs4147359).Except for chromosome 10 SNP (rs4147359), other SNPs have corresponding genes.
According to a previous study, the chromosome 12 SNP (rs3184504) was in the SH2B3 (SH2B adaptor protein 3) gene and is associated with autoimmune disease [38].Multiple studies indicated that SH2B3 was related to the occurrence of autoimmune Hepatitis [39][40][41][42].In addition, recent studies have shown that SH2B3 expressed in lymphocytes might with the risk of mid/long-term clinical relapse after being treated with infliximab in those patients with CD [43].Although how SH2B3 mediates autoimmune disease remains unclear, a study provides us with new insights.Microbiome could exert physiological functions via the SH2B3 gene [44], and gut microbiota also exerts a significant influence on both PSC and UC [45,46].And rs9858213 is in the ring finger protein 123 (RNF123) gene, located in chromosome 3.The protein encoded by this gene displays E3 ubiquitin ligase activity toward the cyclin-dependent kinase inhibitor 1B which is also known as p27 or KIP1, so the research on this gene is mainly focused on tumors now [47,48].A report indicated that p21 expression was higher in IBD cases [49].Unfortunately, no studies have been reported that the relationship between rs9858213 and PSC.
T cells play an important role in both PSC and UC.Many studies focus on T-cell immunotherapy [50][51][52][53].C-type lectin domain containing 16A (CLEC16A) gene which, has been proven associate with multiple immunemediated diseases, which may through T cells to induce pathogenicity [54].This connection validates our results from an immunological perspective.
For rs10909839, this SNP is located in the tetratricopeptide repeat domain 34 (TTC34) gene.TTC34 gene a link with systemic lupus erythematosus was reported by some studies [55,56].Unfortunately, there is limited research on this gene.Therefore, how TTC34 the immune system remains unknown.
We also acknowledge some of the limitations of this study.First, due to data availability, the GWAS data of UC and CD we used were from a European population, while the data of PSC was from a mixed population.In the future, more populations should be included.Second, only 18 SNPs meet the conditions to become IVs.Even if removing linkage disequilibrium, detecting pleiotropy, leave-one-out sensitivity analysis, heterogeneity analysis, and horizontal pleiotropy analysis have been conducted, we cannot guarantee that each SNP site meets the condition that instrumental variables can affect outcomes only through exposure.Some influence of unknown possible confounders inevitably affects our results.We obtained those results by analyzing data from public databases, but the databases didn't provide clinical data.Therefore, experimental or other studies should be conducted to our results.Despite these limitations, our results may inspire possible mechanism analyses and the relationship between PSC and IBD, in the future.

Conclusions
Our study corroborated a causal association between genetically predicted PSC and UC but not for PSC and CD.Then, we identification of shared SNPs for PSC and UC, including rs3184504, rs9858213, rs725613, rs10909839, and rs4147359.

Fig. 1
Fig. 1 Study design and workflow in the present study

Fig. 2
Fig. 2 Scatter plots for MR analyses of the causal effect of PSC on UC and CD.A PSC against UC.B PSC against CD

Table 1
Detailed information on association studies in our analysisUC Ulcerative colitis, CD Crohn's disease, PSC Primary sclerosing cholangitis, SNPs Single nucleotide polymorphisms

Table 2
Detailed information for the SNPs in MR analysis SNPs Single nucleotide polymorphisms, CHR Chromosome, EA EFFECT allele

Table 3
Association of genetically predicted PSC with risk of UC and CD SNPs Single nucleotide polymorphisms, OR Odds ratio, CI Confidence interval

Table 4
The shared SNPs for PSC and UC SNPs Single nucleotide polymorphisms, CHR Chromosome, BP Base pair position, EA Effect allele, OA Other alleles